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ABSTRACT 

We have assembled a sample of 5 X-ray-absorbed and submm-luminous type 1 QSOs 
at 2; 2 which are simultaneously growing their central black holes through accretion 

and forming stars copiously. We present here the analysis of their rest-frame UV to 
submm Spectral Energy Distributions (SEDs), including new Herschel data. Both 
AGN (direct and reprocessed) and Star Formation (SF) emission are needed to model 
their SEDs. From the SEDs and their UV-optical spectra we have estimated the masses 
of their black holes Mbh ^ 10® —10^° ATq, their intrinsic AGN bolometric luminosities 
Lbol ~ (0.8 — 20) X 1O^^L0, Eddington ratios Lbol/LE dd 0.1 —1.1 and bolometric 
corrections Lbol/Lx,2-10 30 — 500. These values are common among optically 

and X-ray-selected type 1 QSOs (except for RX J1249), except for the bolometric 
corrections, which are higher. These objects show very high far-infrared luminosities 
Lfir ^ {2 - 8)xlO^®M0 and Star Formation Rates SFR^ lOOOM0/y. From their 
Lfir and the shape of their FIR-submm emission we have estimated star-forming 
dust masses of Mdust ~ 10® M0. We have found evidence of a tentative correlation 
between the gas column densities of the ionized absorbers detected in X-ray (N//.^^) 
and SFR. Our computed black hole masses are amongst the most massive known. 

Key words: quasars - galaxies: evolution - galaxies: formation - Star formation - 
galaxies: high-redshift - galaxies: starburst. 


1 INTRODUCTION 

In the last two decades, it has become clear that most lo¬ 
cal spheroidal galaxy components (elliptical galaxies and 
the bulges of spiral galaxies) contain a super massive black 
hole (SMBH) in their centres. The proportionality between 
black hole (BH) and spheroid mass suggests a direct link 
between the growth of the black hole as an Active Galactic 
Nucleus (AGN) and the stellar mass of th e spheroid (e.g. 
iMarconi et ahlEopl . iKormendv fc Hdl2013| ~). Identifying the 
main mechanisms for formation and evolution of galaxies, 
and their interrelation to that of the growth of their central 
black holes is a major issue in Astrophysics and Cosmology. 

In the last decade, deep surveys at submillimetre 
(submm) and millimetre (mm) wavelengths have identified a 
high-redshift population of massive dusty galaxies that are 
undergoing extreme starbursts. Since prodigious star for¬ 
mation (SF) is often obscured by dust, these galaxies are 
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luminous in the mm though far-infrared (FIR) wavebands 
where the starlight absorbed by dust grains is re-emitted. 
Likewise, the hotter dust heated in the circumnuclear envi¬ 
ronment of an AGN will emit at mid-infrared (MIR) wave¬ 
lengths. Therefore, submm and MIR observations can be 
combined to study activity in galaxies due to dust-obscured 
starbursts and AGN. The launch of the Herschel Space Ob¬ 
servatory data l|Pilbratt et al.lBoid ') allow obtaining more 
accurate and deeper measurements in the FIR, permitting 
better determinations of the luminosity due to SF {Lsf) and 
hence the SFR. In recent years there have been many studies 
in which the SFR was compared with the growth of the cen- 
tral BH ( Rovilos et al.|[20lAlPage et al.l[2012l . iRosario et ^ 
I2OI2I . iLutz et 20081 ') with conflicting results about the re¬ 
lationship between AGN luminosity and Lsf • 

In this paper, we have stu died a sample o f X- 
ray-obscured QSOs , described by lPage_et^) J2OO4II and 


Stevens et al.l (1 200^1. an d st udied by IStevens et al.l (120041 . 


2010li . Page et al.l I 201lli and iGarrera et all ( 20 111 ') (20 un 


absorbed objects), at z~l-3 when most of the SF and 
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BH g rowth are occurring in the Universe. IStevens et al.l 
ll2005l i found six detections at > 5cr significance at 850/im 
(SCUBA). These QSOs have strong submm emission, much 
higher than typically found in QSOs at similar redshifts 
and luminosities. However for one of them (RX J110431.75 
+ 355208.5) a synchrotron origin for the detected 850/im 
emission cannot be ruled out, so we have not followed up 
RX J110431.75 in this work, and our sample consists of the 
other five sources. They were specifically targeted with Her- 
schel Space Observatory data OT 1 (PI: F.J. Carrera). 

The ultraviolet (UV) and X-ray spectra of our QSOs 
show evidence for strong ionized w inds which produce the 
X-ray obscuration (IPage et al.ll201l]l . Piecing all these clues 
together, we inferred that the host galaxies of these QSOs 
are undergoing strong SF, while the central SMBH are also 
growing through accretion. In pri nciple the ionized winds are 
strong enough to quench the SF JPage et al.ll201ll ') so, given 
that the QSOs have powerful submm emission, these objects 
are then caught at a special time, perhaps emerging from a 
strongly obscured accretion state in an evolutionary stage 
which might last abo ut 10-15 per cent of the QSO lifetime 
dHopkins et al.l [20081 ) . which tallies nicely with X-ray ab¬ 
sorbed broad line AGN bein g about 15 per cent of the X-ra: 


broad line population (e.e. IPage. Mittaz. &: Carreral 


Mateos et al 


lo puiat i 

1 boiOl: 


ICorral et al.l 20 111 : Scott et al.l 


A.-ray 

200d 


20111 : 


Scott. Stewart, fc Mateo? 2012l l. Alternatively, the ionized 


winds are de-coupled from the SF, but then why only X- 
ray-absorbed AGN are submm-luminous remains to be ex¬ 
plained. Some geometrical effect might be invoked, but it is 
difficult to see how to reconcile the dramatically different 
scales at which both processes are expected to happen. 

Here, we endeavour to get the physical properties of the 
central QSOs (luminosities, BH masses, Eddington ratios, 
etc.) and their host galaxies (SFR, Moust, Mqas, etc.) 
and their mutual relationships (or lack thereof). In addition, 
we will try to fathom their place in AGN-host galaxy co¬ 
evolution models. 

The paper is organized as follows. In Section 2, we 
present our sample and summarize the data used. In Section 
3, we present the Spectral Energy Distributions (SEDs) of 
all our objects, explain how we have made fits to several 
representative templates and obtained results from them. 
Moreover, we calculate the Black Hole properties from UV- 
Optical spectra and we study the time scales of the evolution 
of the AGN and the host galaxy. In section 4, we discuss all 
results and compare them with those found in other samples. 
Finally, we summarize our findings in Section 5. 

We have assumed throughout this paper a Hubble con¬ 
stant Ho = 70 km/s/Mpc, and density parameters flm = 0.3 
and Ha = 0.7. 


2 DATA 

We present in Table [T] Herschel Space Observatory data 
for the fields around our QSOs (RX J0057, RX J0941, 
RX J1218, RX J1249, RX J1633). These data include two 
photometric bands (100/rm and lOO/rm) from the PACS 
instrument dPoglitsch et al.l l2010l j and three photometric 
SPIRE bands (250^m, 350^m and 500^m); iGriffin et al.l 
ll2010lf , see Table [2l The PAGS observations were per¬ 
formed in two scan orientations (70°and 110°) for each 


band and field. These images were combined using the stan¬ 
dard pipeline of Herschel Interactive Processing Environ¬ 
ment (HIPE v8 and vlO). The SPIRE were undertaken in 
smallmap mode. We used directly Level 2 data obtained 
from the Herschel Archive and implem ented a script to re¬ 
move the background from the data dUKIRT Newsletter .1 
I 2 OI 0 II . We also present here new optical imaging data (R fil¬ 
ter) from the William Herschel Telescope for RX J1249. The 
observations were taken on 20th May of 2010. The reduction 
and calibr ation were carried o ut using standard iraf proce¬ 
dures fsee lCarrera et ah I [m3 for details). 

In addition, part of the data used in this paper have 
been used or obtained in other studies: X-ray spectroscopic 
obser vations from XMM-New ton and optical-UV spectra 
from iPage et al.| ijgOOlj^ 201j|). Spitzer and SGUBA photo¬ 
metric data from IStevens et M] (I2004l . l^l0ll a nd optical and 
infrar ed photometric data for RX J0941 from ICarrera et ^ 
(l201iri . Finally, we have also used public data from differ¬ 
ent catalogues (SDSS, 2MASS, SUPERCOSMOS, WISE, 
GALEX and OM-SUSS), see Tabled 

We have used Sloan PSF magnitudes, converting them 
to fluxes directly using tuab ~ —2.5 log(Fl,) — 48.60 with 
fvQ = 3.63 X 10® Jy. For u and « we have converted previous 
to AB magnitude using Q: 


uab = usDSS — 0.04 


zab = zsDSS + 0.02 


( 1 ) 


For non-AB magnitudes we have used : 

= F,,o X (2) 

where FL,o are given in column 3 in Table d for each 
band. 

A conservative 5 per cent error was added 
in quadrature to the Herschel (extracted from 
HIPE documentation), SPSS and SUPERCOSMOS 
llHamblv. Irwin, fc MacGillivravl |2001^ catalogued flux 
errors to account for the uncertainties in the zero-points. 
In addition, we have also add ed in quadrature a 2 per cent 
error for the 2MASS data JCohen. Wheaton, fc MegeathI 
I 2 OO 3 II . Regarding data from WISE, we have added a 1.5 per 
cent uncertainty to the catalogued flux errors in all bands 
to account for the overall systematic uncertainty from the 
Vega spectrum in the flux zero-points. Additionally we have 
added a 10 pe r cent uncertaintyto the 12 m and 22 m fluxes 
ilWright et al.l l2010l j to account the existing discrepancy 
between the red and blue calibrators used for the conversion 
from magnitudes to Janskys. 


3 RESULTS 

We have constructed the spectral energy distributions of 
all our objects, in i/Li, versus rest-frame wavelength in pm 
(Figure [T|). We also show the observed SED (from NED) of 
Mrk 231: an ultraluminous infrared galaxy (ULIRG), with- 
o ut any renormal i zation . The average type 1 QSO template 
of iRichards et all 1I2OO6I ) is shown too. They constructed the 
SED from 259 quasars with SDSS and Spitzer photometry, 
supplemented by near-IR, GALEX, VLA and ROSAT data, 
where available. We have rescaled this template to be close 


^ From http://www.sdss3.org/dr8/algoritlims/fluxcal.php 
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Table 1. Summary of the Herschel observations on our fields. 


Object 

RA 

DEC 

z 

OBSID PACS 

Date Obs. 

OBSID SPIRE 

Date Obs. 

RX J0057 

00 : 57 : 34.94 

-27 : 28 : 28.0 

2.19 

1342225345-6 

2011-07-23 

1342234705 

2011-12-18 

RX J0941 

09 : 41 : 44.61 

-1-38 : 54 : 39.1 

1.82 

1342232387-8 

2011-11-17 

1342246614 

2012-06-03 

RX J1218 

12 : 18 : 04.54 

-1-47 : 08 : 51.0 

1.74 

1342233430-1 

2011-12-02 

1342222665 

2011-06-15 

RX J1249 

12 : 49 : 13.85 

-05 : 59 : 19.4 

2.21 

1342235124-5 

2011-12-24 

1342224978 

2011-07-31 

RX J1633 

16 : 33 : 08.59 

-1-57 : 02 : 54.8 

2.80 

1342223963-4 

2011-07-11 

1342219634 

2011-04-26 


Table 2. Photometric data used in this work. All the new data presented in this paper take into account systematic or calibration errors. 
These were added in quadrature to the catalogue errors. 


Filter / Band 

Wavelength 

(pm) 

(Jy) 

RX J0057 
(mjy) 

RX J0941 
(mJy) 

RX J1218 
(mJy) 

RX J1249 

(mJy) 

RX J1633 
(mJy) 

B“ 

0.44 

4270 

0.125 ±0.007 

- 

- 

- 

- 

V* 

0.55 

3670 

- 

- 

- 

0.76 ± 0.04 

- 

R'’ 

0.70 

2840 

- 

0.040 ± 0.003 

- 

- 

- 

Rf 

0.70 

2840 

- 

- 

- 

0.972 ±0.05 

- 

1“ 

0.90 

2250 

0.152 ±0.009 

- 

- 

- 

- 

i'’ 

0.78 

- 

- 

0.072 ± 0.006 

- 

- 

- 


0.3551 

- 

- 

0.0078 ± 0.0010 

0.0130 ±0.0012 

0.391 ±0.02 

0.011 ±0.0014 


0.4686 

- 

- 

0.0175 ± 0.0010 

0.0268 ± 0.0015 

0.715 ±0.04 

0.0365 ± 0.0019 


0.6165 

- 

- 

0.038 ± 0.007 

0.035 ± 0.002 

1.02 ± 0.05 

0.0448 ± 0.0008 

i/c 

0.7481 

- 

- 

0.071 ± 0.004 

0.043 ± 0.003 

1.13 ±0.06 

0.045 ± 0.002 


0.8931 

- 

- 

0.101 ± 0.007 

0.043 ± 0.004 

1.37 ±0.08 

0.052 ± 0.004 

J*’ 

1.25 

1650 

- 

0.14 ±0.02 

- 

- 

- 


2.2 

673 

- 

0.175 ± 0.018 

- 

- 

- 

J'^ 

1.235 

1594 

0.18 ±0.05 

- 

- 

1.55 ±0.07 

- 

H4 

1.662 

1024 

0.36 ±0.07 

- 

- 

1.96 ±0.07 

- 

Rd 

2.159 

667.7 

0.30 ±0.07 

- 

- 

2.54 ±0.11 

- 

WISE19 

3.5 

309.540 

0.306 ±0.014 

0.311 ± 0.012 

0.338 ±0.012 

2.06 ± 0.06 

0.086 ± 0.004 

WISE29 

4.60 

171.787 

0.55 ±0.02 

0.68 ± 0.02 

0.328 ± 0.016 

3.75 ±0.10 

0.101 ± 0.006 

WISE39 

11.56 

31.674 

2.43 ±0.3 

2.4 ±0.3 

1.06 ± 0.16 

17 ±2 

0.23 ± 0.05 

WISE49 

22.09 

8.363 

4.4 ± 1.0 

8.4 ± 1.3 

4.4 ± 1.0 

34 ±4 

0.8 ± 0.3 

IRAC Ch2 

4.5 

- 

0.476 ±0.015 

0.713 ± 0.019 

0.149 ± 0.009 

3.67 ±0.04 

0.083 ± 0.006 

IRAC Ch4 

8 

- 

1.53 ±0.03 

2.13 ± 0.04 

0.363 ±0.015 

12.76 ±0.09 

0.153 ± 0.010 

MIPS 

24 

- 

4.11 ±0.11 

5.25 ± 0.18 

3.28 ±0.15 

26.31 ±0.13 

0.71 ± 0.04 

PACSIOO 

100 

- 

18 ±8 

21 ±7 

11 ±5 

89± 10 

8 ± 5 

PACS160 

160 

- 

31 ± 10 

39 ± 10 

34± 11 

98± 10 

9 ± 6 

SPIRE250 

250 

- 

60 ±7 

89 ±5 

84 ±5 

111 ± 6 

14 ±3 

SPIRE350 

350 

- 

44 ±4 

90 ±5 

71 ±5 

72 ±5 

11 ±3 

SPIRE500 

500 

- 

25 ±4 

68 ±5 

47 ±5 

32 ±4 

11 ±4 

SCUBA450 

450 

- 

32 ± 10 

45 ± 10 

- 

35 ±8 

- 

SCUBA850 

850 

- 

7.8 ± 1.0 

12.4 ± 1.0 

- 

11.0 ±0.7 

6.9 ± 0.7 


^ SUPERCOSMOS. 

^ ICarrera et al. MB- 

SDSS DR7. 

2MASS. 

^ William Herschel Telescope. 
* NED. 

9 WISE All-Sky. 


to the observed SED of our brightest object (RX J1249), 
whose UV-FIR shape it closely matches. 

All objects clearly show at least two components: UV- 
NIR contribution attributable to direct accretion disk emis¬ 
sion (as expected from their type 1 nature), intrinsically ab¬ 
sorbed in the cases of RX J0941 and RX J1218; and a repro¬ 
cessed thermal component in the MIR region from warm op¬ 
tically thick dust further away from the nucleus (the torus). 
In all of them we also observe an additional FIR/submm 
component associated to cooler dust, heated by star for¬ 


mation (SF). We thus confirm the presence of strong FIR 
emission due to SF in these objects, at the ULIRG/HLIRG 
level (compared to e.g. Mrk 231). 

3.1 Fit models 

Our next goal is to make fits to our data with different tem¬ 
plates to extract quantitative information about our objects. 
We use relatively simple empirical and theoretical templates, 
aiming at reproducing the general shape of the SEDs of our 
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Figure 1. SEDs of our five QSOs, compared with a standard QSO template llRichards et al.ll2006l . dashed line, normalized to our most 
luminous object at '--'2000 A) and the observed SED of Mrk 231 (a prototype ULIRG containing an AGN). 


objects. We do not attempt to extract detailed physical in¬ 
formation about our objects from those templates, since our 
data do not warrant such undertaking, and the underlying 
physics is likely to be more complex than that considered in 
the models. We hav e selected the Sherpa: Cl AO’s modelinn 
& fitting package l|Freeman. Doe^ fc Siemig inows B l20nill 
module for the Python platform IIDoe et al.l 2007l l to per¬ 
form the fitting. SEDs were fitted in vLi, versus rest-frame 
A. We have fitted only rest-frame wavelengths longer than 
1218 A, to avoid issues with the Ly-a forest. We have cor¬ 
rected for extinction in our Galaxy using the A(X)/A(V) 
values from Table 21.6 in ICoxI (l200ol i in the observer frame, 
assuming Rv = 3.1 and A{V)/N{H) ~ 5.3 x 10~^^ cm“^, 
extrapolating for wavelengths longer than 250 p,ra assuming 
A{X) oc A-A 

We have modelled the distinct constituents that can 
be appreciated in the SEDs of our objects (see above and 
Figure [T]) with three different components: 

• An AGN accretion disk component: This template mod¬ 
els the direct emission from the accretion disk of the AGN 
(sub-parsec scales). We use the pu r e disk newAGN4 tem¬ 
plate from iRowan-Robinson et al.l (l2008l l. affected by in¬ 
trinsic extinction (see below). The only free parameter is 
the normalization of the template. We have normalized the 
template to its integral in the 0.12-100 pm range (standard 
limits used to estimate the accretion disk luminosity). From 
the fits we directly obtain the value of the disk luminosity 
in that range, Ldisk- 

• A torus component: This template models the re¬ 
emission from the warm and hot dust (o n tens of parsecs 
scales, beyond the sublimation radius, see I Antonuccilll993l l 
that is warmed by the accretion disk emission. We have used 
both an empirical template from iRowan-Robinson et al.l 
(l2008li (dusttor, based on an average quasar spectrum) and 


three dus ty clumpy torus models from lNenkova et al.l (l2008l ) 
found by iRoseboom et al.) (l2013h to represent the average 
properties of t ype 1 QSOs (to r usl, t orus2 and torus3). The 
templates from lNenkova et al.l (l2008ll have Nq ~ 5 — 15 dusty 
clouds along radial equatorial rays. The angular distribution 
of the clouds must have a soft edge, and the radial distribu¬ 
tion decreases as 1/r or 1/r^. 

The difference of the latter model with the others is fun¬ 
damentally that it has a higher inclination of the torus with 
respect to the line of sight (20 deg versus 0 deg), a larger 
number of absorbing clouds (by a factor of three) and a dif¬ 
ferent distribution of the clouds (constant versus declining 
with distance as 1/r). 

Similarly, we have normalized them to their integral in the 
1-300 pm range (again, standard limits used for the torus 
luminosity), so the only free parameter is Ltorus- 


• A SF component: we h ave used a subset of the 
ISiebenm orgen fc Kriieel l (l2007ll sp herical smooth models, 
found by Svmeonidis et alT i 20131 to encompass the ob¬ 
served S EDs of star-forming galaxies at least up to 2 ~ 2. 
The full ISiebenmorgen fc Kriigell (I 2 OOIII models have five 
free parameters to obtain their 7 000 templates. The sub - 
set of templates recommended bv I Svmeonidis et al.l (l2013fl 
(around 2000 templates) are divided into 11 sub-grids, de¬ 
pending on the maximum radius and different star pop¬ 
ulations. For each source, we found the best-fit template 
for each of these 11 sub-grids. We did not attempt to ex¬ 
tract detailed physical information from the particular best- 
fit templates, since instead we were looking for physically- 
motivated templates that reproduced the spectral shape of 
the data. Therefore, the only parameter that we obtained 
from these fits, apart from the best-fit templates, was the 
integrated FIR luminosity (40-500/im) Lfir, since all tem¬ 
plates were normalized in that spectral range. In addition, 
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we have found the optically-thin greybody (GB) model (with 
free temperature T and slope /3) that best matches each 
best-fit SF templates at wavelengths longer than 40 /rm rest- 
frame, providing two more “parameters” associated to each 
fit: T and 

The direct AGN accretion disk component was affected 
by intrinsic extinction modelled in terms of the hydrogen 
column density Nh, so that absorption oc 
where 


r(A) = 


logj, 10 Av A(A) 
2.5 A{V) 


(3) 


We have used Av/Nh = 0.76 x 10“^^ and the 
A{X) /A{V) values from the [Gordon et ah (l2003l i FI SMC re¬ 
sults for A < 8100 A, and from ICoxI ll2000ll for longer wave¬ 
lengths, both parametrizations merging smoothly at that 
wavelength. 

In summary, we have 4 free parameters in the fits: the 
three luminosities (Ldisk, Ltorus and Lfir) and the in¬ 
trinsic extinction (Nh)- In total, we have 44 (4 torus models 
xll SF models) possible combinations of all possible com¬ 
ponents and models. 

We have minimized the usual statistics to choose 
between different models. Attending to the fact that there 
are many “hidden” parameters in our modelling (the dif- 
ferent torus templa tes, and the many parameters i n the 
iNenkova et al.ll2008l and lSiebenmorgen fc Krugelll200^ mod¬ 
els) and to the bad reduced values of our best hts 
(2 — 135), we have neither attempted to use the parame¬ 
ters of the absolute best-htting template as the best value 
of the parameters, nor used Ax^ to estimate the uncertain¬ 
ties in those. Instead, we have taken into account both the 
dispersion of the parameter values among models and the 
relative errors on the photometric points. 

We have inspected the distribution of x^ values for each 
source, finding that there were “families” of models that 
produced quite different values of that statistics. For each 
source (see below) we have chosen the family of models with 
the lowest x^- Table [3] shows the average values of the lumi¬ 
nosities and Nh values among each best-fit family for each 
source and the minimum value of x^ for each best-fit family. 
We have estimated the uncertainties in those values from 
the standard deviation. 

Since have only included among those “families” one 
model for the direct AGN emission and only one or two 
models for the reprocessed AGN emission, the above dis¬ 
persion gives rise to a unreasonably low uncertainty in the 
corresponding luminosities. Hence, we have used the aver¬ 
age relative error on the relevant photometric points to es¬ 
timate an additional relative uncertainty in all luminosities. 
We have calculated the average relative error in the rest- 
frame ranges 1218-10000 A, 3-30 /rm and > 100/rm for the 
direct AGN, reprocessed AGN and SF components, respec¬ 
tively. Those average relative values were multiplied by the 
corresponding average luminosities, finally adding them in 
quadrature to the above standard deviation. 

We believe that these values and their uncertainties are 
fair estimates of the luminosities of each component and 

^ We obtained the numerical values from 
http: //www. stsci. edu/~kgordoii/Dust/Extinctioii/MC_Ext/nic_ave, 


of the effects of the photometric errors and our lack of an 
accurate physical model for what is really happening in each 
source. 

The shap e of most of our SED s in Figure [T] is very 
similar to the [Richards et al.l ([2006[ f (R06) QSO template 
(UV-MIR region), but with an excess in the region of star 
formation (far-IR /submm region). So we have repeated 
the fits replacing the obscured disk and torus components, 
with intrinsic extinct ion, type 1 QSO SED template from 
[Richards et al.[ l|2006h . keeping the SF component (a total 
of 11 sub-grids per fit). We have obtained x^ values and FIR 
luminosities which are quite similar and compatible with the 
previous three-component ones. From these results we can 
draw two main ideas: 

• The shape of our objects is well modelled (except per - 
haps for RX J1633, see below) by the [Richards et al.[ (|2006[ f 
template (UV-MIR region) and SF component, with some 
intrinsic extinction in two cases: they do not appear to be 
unusual in the UV-MIR region. Fits with separate disk and 
torus components are similarly good. In what follows we 
quote the latter results because they allow a more straight¬ 
forward estimation of the direct and reprocessed AGN emis¬ 
sion, and because they are clearly better for RX J1633. 

• There is an excess in the far-IR /submm region that 
does not look like emission from the torus, indeed it is fitted 
very well with SF templates. We obtain quite similar FIR 
lu minosities using accret ion disk and torus models together 
or [Richards et aP (|2nn6ll template. 

On the other hand, all of our objects (except RX J1633) 
prefer the torus models dusttor, torusl and torus2 which 
have similar shape. However RX J1633 clearly prefers the 
toru s3 model. This t orus model is one of the most extreme 
from [Nenkova et al.[ ([2008l f with respect to the cold gas emis¬ 
sion: it has a peak at higher wavelengths, but it is not enough 
to model our excess in far-IR /submm region. Summariz¬ 
ing, this (together with the discussion at the beginning of 
the Section 13.511 reinforces our idea that the excess far-IR 
/submm region corresponds to SF emission and not to the 
reprocessed AGN emission. 

In the following, we have used the full three compo¬ 
nent fits for all QSOs because they treat the direct and in¬ 
direct AGN components separately, allowing a comparison 
between them. 


3.2 Luminosities from SED fits 

We show in Table [3] the average value of the 8-1000 /rm lu¬ 
minosity Liu, calculated by integrating the torus and SF 
templates in that range and adding both contributions. Its 
uncertainty has been calculated in a similar way to those 
of the best-fit luminosities (see above), using both the stan¬ 
dard deviation among the best-fit families and the relevant 
photometric errors (using this time the 8-1000 /rm range), 
adding them in quadrature. 

We define Lbol as the integrated 100 keV-to-100 /rm 
luminosity from the AGN. We estimate this quantity us¬ 
ing the observed absorption-co rrected 2-10 keV (obtained 
from XMM-Newton spectra in [Page et al.[ lioill j and 0.12- 
100 /rm luminosities of our objects. For wavelengths longer 
.extthainiI.12 /rm we have integrated the fitted disk template. At 
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rxj0057 SED 



1 10 

Resl-frame wavelength (n m) 


100 



1 10 100 
Rest-frame wavelength (^l m) 



1 10 100 
Rest-frame wavelengthfji m) 



1 10 100 
Rest-frame wavelength (^l m) 


rxj1218 SED 



rxj1218 SED 



Figure 2. Best-fit SEDs to each source (one row per source) in two columns: on the left separate direct AGN (green dashed line), 
reprocessed AGN (torus, orange dashed line) and SF (blue dashed line) components; on the right, total AGN R06 (green dashed line) 
and SF (blue dashed line) components. We also indicate for each source which is the best-fit torus model. First row: RX J0057: Torusl. 
Second row: RX J0941: Dusttor. Third row: RX J1218: Torus2. Fourth row: RX J1249: Torusl. Fifth row: RX J1633: Torus3. 


wavelengths shorter than 0.0035 fim (0.35 keV) we have ap¬ 
proximated the source continuum emission by a power-law 
Li, oc normalized to the observed absorption-corrected 
2-10 keV luminosities ijPage et al ] l201lh . Between those two 
wavelengths we have interpolated linearly in log-log (equiva¬ 


lent to assuming a power-law shape). Most of the final bolo- 
metric luminosity comes from the disk component (50-70 
percent, except for RX J1633, for which this is 34 percent), 
so the exact parametrization of the higher energies does not 
have a strong influence on our results. We have calculated 
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Figure 2 - continued 


Table 3. A verage results and their uncertainties of fit parameters for each QSO, except for the X-ray AGN luminosity Lx, AGN which 
comes from l^ge et aul201ll . represents the minimum value for each best-fit family and N represents the number of photometric 
points that we have for each object. 


Object 

Lx,AGN 

(10“L©) 

Nh 

(10^° cm-2) 

Ldisk 

(IQiiLo) 

LtORUS 

(IQiiL©) 

Lfir 

(iQiiLe) 

Lir 

(IQiiL©) 

T 

(K) 

/3 

XVN 

RX J0057 

2.93 

0.01 ±0.05 

125 ± 17 

190 ± 30 

49 ± 10 

170 ± 40 

38 ±4 

1.1 ±0.2 

159 / 20 

RX J0941 

0.93 

0.77 ±0.11 

100 ± 10 

150 ± 20 

78± 11 

180 ± 30 

36 ±3 

1.2 ±0.2 

347 / 23 

RX J1218 

2.33 

0.39 ±0.08 

42 ±7 

25 ± 5 

63 ±6 

90± 19 

36.0 ± 1.7 

1.2 ±0.2 

414 / 17 

RX J1249 

3.69 

0.14 ±0.03 

1230 ± 90 

1160 ± 100 

70± 15 

780 ± 90 

39 ±4 

1.2 ±0.2 

3109 / 23 

RX J1633 

5.85 

0.003 ± 0.002 

65 ±6 

70 ±30 

22 ± 5 

90 ±30 

36 ±2 

1.00 ±0.11 

29 / 17 

Table 4. Additional physical quantities derived from the average values and their uncertainties of the best-fit parameters. 


Object 

Lbol 

(IQiiL©) 

CF 

SFR 

(M01/-1) 

MdUST 
(10® Mo) 

^EDD 

M 

ksOL 

RX J0057 

250 ± 30 

0.78 ±0.15 

840 ±170 

2.5 ±0.8 

0.09l®;°5 

17 ±2 

85 ± 10 

RX J0941 

169 ± 13 

0.87 ±0.14 

1350 ± 190 

4.1 ± 1.3 

0.09t°;°| 

11.5 ±0.9 

182 ± 14 

RX J1218 

78 ±7 

0.32 ±0.07 

1090 ± 100 

4.6 ± 1.6 

0.l3l°;°2 

5.3 ±0.5 

33 ±3 

RX J1249 

1890 ± 90 

0.61 ±0.06 

1200 ± 300 

3.8 ±0.9 

d d 
+ 1 

CD 

128 ±6 

512 ±8 

RX J1633 

186 ± 15 

0.40 ±0.16 

380 ± 90 

2.0 ±0.7 

1 -,+0.6 

12.6 ± 1.0 

32 ±3 
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Table 5. Black hole masses from fits to the CIV line at 
1548A(except for the row marked with *, which comes from a 
fit to Mgll at 2798A), FWHM corresponding with each line, 
the monochromatic luminosity corresponding with each line and 
the Eddington luminosity. 


Object 

FWHM 

(A) 

APa 

(lO'^'^ erg/s) 

log{MBH) 

(log(M0) 

^OgiLEDo) 

(log(T0)) 

RX J0057 

146+3 

419+ 18 

9.94+0.36 

14.45+0.36 

RX J0941 

113+24 

322 ±29 

9.77+0.40 

14.27+0.40 

RX J0941* 

106+7 

202 ±18 

9.22+0.55 

13.72+0.55 

RX J1218 

77+23 

150 ±28 

9.28+0.45 

13.79+0.45 

RX J1249 

86+26 

3731+392 

9.99+0.45 

14.49+0.45 

RX J1633 

53.9+1.2 

178 ±6 

8.73+0.36 

13.24+0.36 


Lbol for each best-fit family, showing in Table |4] their av¬ 
erage and its uncertainty, coming again from the dispersion 
among the best-fit families and the average relative photo¬ 
metric errors (using the same range as the one for the direct 
AGN component), adding them in quadrature. 


3.3 Black Hole masses and related quantities from 
optical spectral fits 

We have estimated the BH masses of our objects from the 
M gll and CIV emis sion lines from the optical-UV spectra 
in IPage et al.l (l201ll i. The latter is usually preferred in the 
literature over the former because it presents lower complex¬ 
ity and wider scatter. Unfortunately, Mgll is only covered 
by our spectra for RX J0941; for this source we have es¬ 
timated the BH mass using both lines, to gauge the dif¬ 
ference between using them. There are a number of dif¬ 
ferent parametrizations of the BH mass as a function of 
the width of those lines and the luminosity in neighbour¬ 
ing w avelengths. We have chosen those used in IShen et al.l 

1I2OIII) : 


iog(^^^2iogr.^^u 


Me J 

+0.50 log 


V1000 km/s ) 
ALa 


+ 6.86 


^10** erg/s^ 
for Mgll llVestergaard fc Osmeillioogli . and 


(4) 


log 


/ Mbh \ 


2 log 


/ FWHM \ 
V1000 km/s ) 


+ 


+0.53 log 


ALa \ 
1044 erg/s J 


+ 6.66 


(5) 


for CIV llVestergaard fc Petersonll2006l i. with rms values of 
0.55 and 0.36 dex, respectively. 

To obtain the FWHM of the Mgll and CIV lines, we 
have followed the technique outlined in Shen et al. (2011), 
which we summarize here. For Mgll we have fitted, over the 
2700-2900 A rest-frame range, one narrow Gaussian com¬ 
ponent with FWHM (km/s) = ma x(AA, 500) w here A A is 
the spectral resolution as given in IPage et al.l (1201 il l and 
one broad component, plus four Fe “humps” at rest-frame 
wavelengths of 2630, 2740, 2886 and 2950 A. The continuum 
has been modelled as the best-fit power-law over the rest- 
frame range 2200-2700 A (since the additional 2900-2090 A 



MA) 


Figure 3. An example of the fits made to estimate the FWHM 
of the CIV emission line. We have modelled this feature for 
RX J0057 with three additive Gaussians (blue, orange and vio¬ 
let dashed lines), a power-law continuum (dotted green line) and 
three absorbing Gaussians (not shown individually). The total 
best-fit model is shown as a red solid line. 


range is not covered in our spectrum of RX J0941). For CIV 
we have fitted three Gaussians (only one was necessary for 
RX J0941) over the range corresponding to the expected 
position of CIV at the redshift of the source ±20000 km/s, 
forcing their central wavelengths to be in the range 1500- 
1600 A. The continuum has been modelled as the best-fit 
power-law over the rest-frame ranges 1445-1465 and 1700- 
1705 A, see Figure [3] 

In addition, the spectra of our sources presented a num¬ 
ber of narrow and broad absorption lines, which we have 
modelled as multiplicative absorbing Gaussians M{\) = 
^1 — Aexp with 0 < A < 1. We needed 3, 

5, 7, 3 and 6 such components for RX J0057, RX J0941, 
RX J1218, RX J1249 and RX J1633 (respectively) around 
CIV, and 2 for RX J0941 around Mgll. 

For RX J0941 we have obtained the 1 — a uncertainty 
on the FWHM directly from that of the a of the single 
broad line fitted to the Mgll and CIV lines, obtained in 
turn from the Ax^ = 1 interval around the best fit value. 
Unfortunately, for the rest of the sources this is not so 
straightforward, since the CIV emission is modelled as the 
sum of three Gaussians, and there is no analytic expres¬ 
sion relating the overall FWHM to the best fit param¬ 
eters of those three Gaussians (nine in total: three each 
of central wavelength, width and normalization). We have 
therefore sampled the parameter space, fixing those nine 
parameters to values around the best-fit ones, calculat¬ 
ing the best-fit to each spectrum in the above range with 
the above constraints, calculating the FWHMi for each of 
those combinations of parameters, and assigning to each a 
probability Pi = (Press et al. 1992). We have 

then estimated the variance on the FWHM as ctfwhm ~ 
{FWHM - {FWHM)) ~ Y.i Pl{FWHMi - {FWHM))^ 
where {FWHM) ~ Y.i Pi FWHMi and P/ = Pi/ Yii Pi (to 
normalize the total probability to unity). For the sampling 
of each parameter we have used a Gaussian centred on the 
best-fit value with a dispersion set to (in decreasing order of 
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choice) either the corresponding value of the covariance ma¬ 
trix (if defined), or to the Ax^ = 1 interval (if defined), or to 
10% of the best-fit value (1% in the case of the central wave¬ 
lengths). The latter percentages were the rough averages of 
the dispersions of the similar parameters when they were 
defined. In any case, we checked that the final result did not 
depend on the particular choice of these dispersions, since 
they only serve to increase the efficiency of the sampling, 
favouring values close to the best-fit, where the highest prob¬ 
abilities should concentrate. We sampled 2000 combinations 
of parameters for each spectrum. We also checked that the 
dispersions did not depend on the exact number of sampling 
points. 

The estimates of the continuum values vLi, in the Eqs. 
m [S] (at rest-frame wavelengths of 1350 and SOOOA for CIV 
and Mgll, respectively) have been obtained from the disk 
component that best fitted the overall SED of each source 
(see Section \m . which should provide a good estimate of 
the intrinsic disk emission, corrected for absorption. We have 
estimated the 1 — a uncertainties on those values from the 
dispersion in the values of the normalization of that compo¬ 
nent Ldisk (see Section [3Tl) . 

Finally, the uncertainties in log{MBH /Mq) have been 
estimated from those of FWHM and of the continuum using 
the standard error propagatio n rules on the equations above 
dBevington fc RobinsonI Il992l l. We have added in quadra¬ 
ture the rms values in the equations above to get our rather 
conservative estimates of the total uncertainties. The BH 
mass estimates and the total uncertainties are given in Ta- 
ble[S] For RX J0941 we get log^MsH /Mq) = 9.8 ±0.4 from 
CIV and log^MsH /Mq) = 9.2 ± 0.6 from Mgll, which are 
within their mutual errors. 

From the mass of the black hole we can calculate the 
Eddington luminosity using l^dingtonl lll913ll : 

Tbdd = 3.2x1O^^^L0 (6) 

Mq 

3.4 AGN-related quantities from SED fits 

Once we have Lbol we can estimate naively the covering 
factor. Due to the anisotropy of the torus radiation we can 
only estimate t he bolometric covering factor or apparent 
cover ing factor (iMor. Netzer. fc Elitzuil l2009l : iHonig et al.l 
l2014l 'l CF as the ratio between the reprocessed emission from 
the torus and the total AGN bolometric luminosity: 

CF=i^tp^^ ( 7 ) 

Lbol 

This gives us a measure of how much of the sky seen 
from the central source is intercepted by the torus. The val¬ 
ues, shown in Tabled) have been calculated from the average 
values of Ltorus and Lbol and their uncertainties, using 
the standard propagation of errors. 

The Eddington ratio would then be Aedd = 
Lbol/Ledd- We show the values of Aedd in Table ID they 
span the range ~ 0.1 — 1.1. The error bars on Lbol are 
symmetric in linear space, while those of log Mbh are large 
and symmetric in log space. The usual propagation of er¬ 
rors rules are only formally valid for small errors and using 
them would result in symmetric error bars for Aedd larger 
than the actual values. This does not make sense, since both 
the bolometric luminosity and the black hole mass are very 


significantly detected. To avoid these difficulties, we have es¬ 
timated the uncertainties on the Eddington ratios sampling 
10000 times the distribution of Lbol and log Mbh using 
Gaussians with the observed values, calculating the corre¬ 
sponding Eddington ratios, and finding the narrowest inter¬ 
val that included 68.3 percent of the sampled values. The 
resulting asymmetric error bars are those shown in Table [D 


3.5 Star-formation-related quantities from SED 
fits 

Our SEDs require a contribution from cold dust, quantified 
by Lfir- We have used the luminosities obtained in the pre¬ 
vious Sections (given in tables 13.21 and [3| to check whether 
that cold dust could be powered by the AGN, following a 
simple energetic argument. We have compared Lbol (the 
total AGN power available) with Ltorus (AGN power al¬ 
ready intercepted by the "warm” reprocessing medium -the 
torus-) and Lfir (the observed power from cold dust): if 
Lbol is significantly larger than the sum of Ltorus and 
Lfir, the AGN would have enough power to generate all 
the observed emission (with some to spare for the, mostly 
unobscured, observed UV-to-optical range). This can also be 
verified visually by looking at the SEDs in Fig. O It is clear 
that the AGN emission would be insufficient in RX J0057, 
RX J0941 and RX J1218. We have already discussed that the 
preferred torus model for RX J1633 includes a larger contri¬ 
bution from colder gas, producing a lower LFIR, so some of 
the putative AGN contribution would be already corrected 
for. As we will see below, RX J1249 is an exceptionally lu¬ 
minous QSO. However, its FIR emission is comparable to 
that of the other members of our sample, which would be 
unexpected it a significant contribution to that range came 
from such a powerful AGN. Concluding, the FIR emission 
observed in our objects likely predominantly comes from SF, 
as obviously it’s not just star-formation. 

The SFR has been calculated from Lfir using the fol¬ 
lowing expression from iKennicuttI (ll998l L 

SFR{MQ/y) = 1.7217 x 10 -^“Lfir(To), (8) 

We obtain values of SFR ~ WOOMq /y (shown in Ta¬ 
ble [3|. We confirm then our qualitative impression from 
Fi&m these objects are forming stars copiously. 

Supernovae and bin ary stars associated with star for¬ 
mation produce X-rays. ISvmeonidis et aP (1201 ill found a 
tight correlation between X-ray luminosity and IR lumi¬ 
nosity from ULIRGs, which can be expressed as Lx,sf = 
(1.9 X 10^®)Lj;r ± (4.15 X 10“®)Lji? (their equation 3). For 
the typical values of our QSOs of Lm ~ 10^®, this cor¬ 
responds to Lx,sf ~ erg/s. Comparing these val¬ 

ues with the AGN X-ray luminosities in the same band 
(Lx-agn ~ 4 — 20 X 10^"* erg/s), we conclude that the 
contribution of the SF to the X-ray luminosity is very small 
and can be neglected. 

In order to estimate the values of the dust mass as¬ 
sociated with star formation, we have used the values of 
Lfir, and the T and /3 obtained in Section 13.11 for each 
of the “good” sub-grids. We have related those three quan- 
ti ties to the dust mass Mdust using Equations 2 and 5 
in lMartfnez-Sansigre et ^ (l2009l l. Their Equation 2 can be 
written as 
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S-Kh dl{z)F^, 1 

c2 il + z){i^'{l + z))PB^{T,v'{l + z)) F 


where F^,/ is the observed monochromatic flux at some ref¬ 
erence observed frequency u' (in their case A' = 1.2 mm), 
where 






dv 


r(/3 + 4)C(/3 + 4) 


exp{hu/kT) — 1 


f h 
\kTJ 


( 10 ) 


( llGradshtevn et al.ll2007ll FI II 792a) is the normalization of 
the greybody function to the full frequency range. Eq. [5] is 
simply the total luminosity of a greybody normalized to the 
observed value of F^i. We have called this quantity L'pjn to 
emphasize the difference between this and the usual defini¬ 
tion of Lfir over the 4 0 to 500 /xm range. _ _ 

Similarly, Eq. 5 in iMartfnez-Sansigre et al.l ll2009ll can 
be written as 


A7dust = 


_ dl{z)F,' _ 

(1 -I- z)k,{i''{1 + z))B^{T, !z'(l -I- z)) ’ 


( 11 ) 


where 


k(iz(1 + z)) = ko X . (12) 

We h ave used kq = 0.04 m^/kg at A' =1.2 mm jBeelen et al.l 
l200(l l. If we now divide Eq.[TT]by Eq.[9l substituting Eg. 1101 
we get 


Mdust = 


c" 1 
Stt/i KoU'~P 


L'p 


POO j i/O-tP 

Jo ^ eyip(hu/kT) — l 


(13) 


where the last fraction is simply a normalization over the 
full frequency range. 

We have instead used 
^2 1 


Mdust = 


Snhuou'-f^ 

J 40 him 


t,3+)3 


(14) 


eyip{hiy / kT) — 1 


using the usual definition of Lfir, which we have obtained 
directly from our SED fits. We show in Table |4] the aver¬ 
ages and standard deviations of the acceptable families of 
models. Our approach involves two approximations: replac¬ 
ing the renormalization over the full interval to that over 
the 40 to 500 fim interval, and replacing the FIR luminosity 
from a greybody with that from a fit to a more sophisticated 
SF model. We believe that we are justified to do both be¬ 
cause, on the one hand, the greybody shape decreases very 
fast outside the usual FIR range and, on the other hand, 
the integrals over that range of the greybody fits to the SF 
templates gave very similar values to those of the actual 
templates. 

We can also estimate the gas mass present in the star¬ 
forming regions of those galaxies from the du st masses, as¬ 
sumin g a gas-to-dust ratio of 54, deduced by iKovacs et al.l 
(l2006ll for z = 1 — 3 SMG (with an uncertainty of about 
20 per cent). Th i s valu e is similar to the one obtained 
by ISeaauist et al.l (l2004h for the central regions of nearby 
submm bright galaxies. For our typical dust mass of ~ 
10 ®Mq this corresponds to a gas mass of Mqas ~ 10^® — 
10“ Mq. 


3.6 Time scales 

Assuming an efficiency in conversion of ac creted mass into 
radiation 77 (we have assumed 10% e.g. iTreister fc UrrvI 


Table 6. Summary of timescales for QSOs and star formation, 
assuming constant mass accretion rates and SFR. r is the black 
hole mass-doubling timescale, Tmax is the time needed to reach 
the maximum local black hole mass (see text), and trf is the 
time needed to reach the corresponding maximum host galaxy 
mass, t is the “look-back time” (see text). 


Object 

T 

'^max 

TSB 

t 


(Gy) 

(Gy) 

(Gy) 

(Gy) 

RX J0057 

Q c+0.3 
U.O_o.4 

'-0.4 

10 ±2 

10.6 

RX J0941 

A c+O.S 

1 q+0.5 
-*-•^-0.4 

6.3 ±0.9 

10.0 

RX J1218 


0 4+0.5 

7.8 ±0.7 

9.9 

RX J1249 

0.08t°°® 

O-OStHe 

7.1 ± 1.8 

10.6 

RX J1633 

0-041^4 

1-551S+ 

22 ± 5 

11.3 


I2OO6II the mass accretion rate M can be related to the bolo- 
metric luminosity Lbol by: 

Lbol = rjMc^ (15) 


The relative speed of galaxy-building through star formation 
and black hole growth through accretion can be estimated 
from SFR/M. Both quantities are included in Tabled with 
errors derived using the standard propagation of errors. 

Assuming a constant accretion rate, we can estimate 
the black hole mass doubling time r as 


Mbh 

M 


(16) 


and, defining the maximum mass of the black hole as 
Mbh, max (which we take as 2 x 10^® Mq, the maximum 
value of the observed black hole masses in t he local Universe, 
NGC 4889, from coluinn 10 in Table 2 of iKormendv fc Hoi 
l2013l . [McConnell et al.ll2012h . the time Tmax needed to reach 
it is 

Mbh ,max Mbh 

'T~max :- v-^0 

M 

Alternatively, we could have used in the numerator of Eg. 1171 
simply Mbh, max, but the masses of the black holes in the 
centres of our objects are already comparable to the maxi¬ 
mum local black hole mass, so Eq. E] is more accurate. We 
have again sampled the distribution of M and log Mbh to 
estimate the uncertainties in these two timescales, as dis¬ 
cussed above. 

Assuming again a maximum local black hole mass, we 
can estimate what would be its corresponding maximum 
hos t galaxy mass, MBULOE .max = 8.53 x lO’^® Mq, using 
the [Marconi fc Hunt! (I2OO3I I relation. From that value and 
the current SFR observed in our objects, assuming again 
constant SFR, we can estimate the time tsb needed to reach 
that maximum mass value as: 


tsb 


Mbulge ,1 
SFR 


(18) 


The values of all these timescales are shown in Table E] along 
with the time between their epochs and now (“look-back 

timh’El)- 


^ From http://www.astro.ucla.edu/wright/CosinoCalc.htinl 
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4 DISCUSSION 


In this Section we will piece together the clues obtained 
above about the nature of our objects. 

From the SED shapes and template fits, we see some 
differences in the properties of the objects in our sample: 
RX J0057 and RX J1249 have similar redshifts z ~ 2.2, 
small optical obscurations, the highest bolometric AGN lu- 
minosities and a prefe rence for the first torus model in 
iRoseboom et al.) (l2013tl . RX J0941 and RX J1218 have 
again similar redshifts 2 ~ 1.8, higher optical obscurations, 
the lowest b olometric AGN luminosi t ies an d also allow for 
the empirical Rowan-Robinson et al.l (l2008l i and the second 
IRoseboom et al.l ( 2013l i torus models. Finally, RX J1633 has 
the highest redshift 2 = 2.8, small optical obscuration, in¬ 
termediate bolometric AGN lumin osity and a strong pre f- 
erence for the third torus model in IRoseboom et ^ (l2013li . 
This results in a higher torus contribution at longer wave¬ 
lengths, intuitively corresponding to a higher probability of 
further reprocessing of the direct emission at larger distances 
and lower temperatures. Additional evidence for a (relative) 
higher inclination of RX J1633 comes from its radio mor¬ 
phology in FIRST images, which shows two diffuse blobs 
roughly symmetrical with respect to the optical position of 
the QSO, in contrast with the rest of our objects, which are 
either not detected (RX J0057 and RX J1218) or pointlike 
(RX J0941 and RX J1249). 

The host galaxies of RX J0941 and RX J1218 seem 
to be growing much faster than their black holes, both in 
absolute terms and compared to the other objects. This is 
because their AGN are less luminous but have similar SFR 
to the rest. The difference in AGN luminosities between the 
objects could be ascribed to the redshift, since those further 
away are also more luminous, as expected in a flux-limited 
sample, but this is belied by the intermediate luminosity of 
the highest redshift object RX J1633, which is also the one 
with the most discordant torus properties. In any case, with 
such a small sample it is impossible to assess the significance 
and implications of the differences observed. 


4.1 AGN properties 


In this section we want to study the main properties of 
the AGN. Taking the sample as a whole, we can compare 
their overall properties to those of type 1 QSO from the 
SPS S DR7 dAb azaiian et al.l l2009l j: bolometric luminosi- 
ties (IShe n et al.r 201ll 'l. Eddington lumino sities and ratios 
jKgllv &: Shen 20131 ') and covering factors (IRoseboom et all 

I2OI3I) . 


As expected from the design of the sample 
dPage. Mittaz. fc Carrer^ |200(]| ~). their X-ray and bolo¬ 
metric luminosities are similar to those of QSOs in their 
redshift ranges (using a range A 2 = 0.4), except for 
RX J1249, which is truly exceptional: there are only 50 
objects br ighter than RX J1249 (<0.05%) in the whole 
sample of I Shen et al.l (l201ll 'l. and it is one of the most 
luminous object even compari ng with the high-b olometric- 
luminosity-selected sample of iTsai et aP (l2014l 'l. Also, as 
discussed above, the shape of their SEDs in the optical-MIR 
region is unexceptional. 

From Tabled three of our QSOs have Eddington ratios 
around 0.1 (RX J0057, RX J0941 and RX J1218), one has 


a value around 0.6 (RX J1249) and the last one is closer to 
unity (RX J1633), although with large error bars, specially 
in the las t two. This range of v alues is compatible with that 
found bv iKellv fc Sheiil (120131 '). who have studied a sample 
of ~ 58000 Type 1 QSO, also finding that Aedd ^0.1 are 
significantly rarer at 2 ~ 4 compared to 2 < 2, while objects 
with Aedd ^0.1 are found similarly in both redshift ranges. 
This is also similar to what we see in our sample, where the 
four objects with 2^2 are compatible with Aedd ~ 0.1, 
while RX J1633 at 2 = 2.8 shows some preference for higher 
values of Aedd- They conclude that Type 1 quasars radiat¬ 
ing near the Eddington limited are extremely rare, suggest¬ 
ing that Type 1 quasars violating the Eddington limit do so 
only for a very brief period of time. RX J1633 could be in 
this situation. _ 

IRoseboom et ahl (1201311 have estimated the covering fac¬ 
tors CF of their WISE-UKIDSS-SDSS (’WUS) quasar sam¬ 
ple (5281 quasars) with 2 < 1.5. They have found an av¬ 
erage value of CF = 0.39, but in addition they also get 
that two-thirds of type 1 quasars have CF in the range 0.25 
to 0.61, roughly as in our sample. Segregating the sources 
according to their bolometric luminosities, our four lower 
luminosity objects are in their 46 < log(I/soL) < 47 bin. 
A quick MonteCarlo sampling of the log-normal fit to their 
CF distribution (calculating values from our sample with 
respect to the mean values of the distribution and hve ran¬ 
dom distribution values with respect to the mean values too) 
shows than an unremarkable 40 per cent of the simulated 
samples have a similar distribution as our objects in this 
bin. In contrast, RX J1249 is in their highest luminosity bin 
(log(Lsoi.) ~ 48), where only three percent of the objects 
have a higher CF than this source: RX J1249 is not only 
one of the most luminous objects known, it also has an ex¬ 
ceptionally high reprocessed-to-bolometric AGN luminosity 
ratio for its luminosity. 

We have calculated the bolometric corrections to the 
X-ray luminosities for our sample dehned as: 

kBOL = ( 19 ) 

-^A',2-10 

finding values ksoL ~ 30 — 500 (see Table 0 for Aedd ~ 
0.1 — 1. Although at face value this lo o ks ve ry different 
from the results of IVasudevan &: FabiarJ (l2009ll , our three 
objects with Aedd ~ 0.1 (RX J0057, RX J0941 and 
RX J1218) are well within the values spanned by the ob¬ 
jects in their sample at similar Eddington ratios ksoL ~ 
30 — 180 (see their Fig. 6). The most discrepant objects 
are RX J1249 and RX J1633. The former presents an ex¬ 
tremely high bolometric correction around 500, more in 
line with its BAL nat ure (iGruoe. Leighlv. fc Komossall2008l : 
iMorabito et al.l[2014^ . On the contrary, RX J1633, our high¬ 
est redshift object, presents a bolometric correction of about 
30 for an Eddington r atio o f about 1. Only one object from 
IVasudevan fc Fabiani (l2009l j has a smaller correction for that 
ratio, although there are a few more with Aedd > 0.1 
and ksoL < 40 in their sample, so it is not exceptional. 
We have also compared our bolometric corrections with 
those of iMarconi et all (l2004l'l . Only RX J0057 is close to 
their relation and inside the dispersion region, RX J1249 
and RX J0941 are above the relation and RX J1218 and 
RX J1633 are below the curve. Again the most discrepant 
sources are RX J1249 and RX J1633. 
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4.2 Star formation properties 

Regarding now the star formation properties of our sam¬ 
ple, we have obtained very high SFRs~ IOOOMq /y: these 
objects are forming stars copiously and they are at the 
ULIRG/HLIRG level with strong FIR emission. 

The distribution of greybody temperatures (T = 36 — 
39 K) of the SF templates of our s ources seems to b e 
closer to those of SMGs (T ~ 35 K. iBeelen et al.ll2006ll . 
than to those of other high redshift QSOs (T ~ 47 K, 
iKovacs et al ][200i): the range of temperatures in the lat¬ 
ter work is T 40 — 60 K while only two of our sources 
have temperatures compatible with the lower bound of the 
interval (RX J0057 and RX J1249, the first subset discussed 
at the beginning of this section). On the other hand, we 
obtain greybody slo pes (/9 = 1 — 1.2) below those of typ¬ 
ical SMG (/3 = 1.5 iBee len et a ll l2006l l and high redshift 
QSOs (/? = 1.6 iKovacs et al.l 2006h . even taking into ac¬ 
count the substantial error bars. However, our temperatures 
and slopes could be affected by the limited wavelength range 
of the greybody fits to the SF templates, since the SF ther¬ 
mal bump in our sources is generally broader than in simple 
greybody models, and hence the fit tends towards low values 
of T (to match the steep decline at the longest wavelengths) 
and low values of /3, to accommodate the extra width at 
the shorter wavelengths. If we just use a simple greybody 
model to pa rametrize the SF, instead of the more elaborate 
templates of ISiebenmorgen fc Kriig^ (l2007l l , we confirm the 
above tendency: both RX J0057 and RX J1249 prefer higher 
temperatures than in the limited-range-fit to the best-fit 
SF template, while the rest show similar or slightly lower 
temperatures. This would go in line with the fact that in 
the former two objects the maximum of the SF emission is 
lower than that of the torus emission, while in RX J0941 
and RX J1218 the opposite happens. Finally, it might well 
happen that a single temperature fit without some modifi¬ 
cation at the high frequency end is not sufficiently accurate 
to allow a precise determination of /3-T values. In any case, 
the exact values of those parameters in particular fits have 
a very limited impact in our inferred luminosities, since we 
derive quantities from averages and dispersions over the best 
group of fits. 


4.3 AGN-SF relationship 

In the context of possible evidence for an influence of the 
central AGN on the evolution of its host galaxy, it is in¬ 
teresting to compare the rate of black hole growth (as 
gauged from e.g. the X-ray luminosity Lx, 2 - 10 , the intrin¬ 
sic bolometric luminosity from AGN Lbol, or the accretion 
rate M) to the rate of galaxy growth (from th e SFR, the 
LIR o r the monochromatic 60/im luminosity). iLutz et 'aD 
ll2008tl studied the correlation between log{uL,^) at 60/rm 
and PAH emission (as proxies of the SF luminosity) and 
log{uLi,) at 5100A(as a proxy for the AGN luminosity) for 
a sample of QSOs at similar 2 ; to ours with Spitzer spec¬ 
troscopic data in the rest-frame mid-infrared. They have 
found that, at high luminosities and 2 , there was a flat¬ 
tening of the relation between SF and AGN luminosity 
that i s observed for lower redshift QSOs. iMullanev et al.l 
ll2012lf used 100 and 160 /rm fluxes from GOODS-Herschel 
finding no evidence of any correlation between the X-ray 


and infrared lu minosities of modera te luminosity AGNs 
at any redshift. iRovilos et al.l ll2012lf found a significant 
(>99%) correlation between Lx, 2-10 and SFR (see Fig¬ 
ure [ 4 ]) in a deep GOODS-XMM-Newton-Herschel sample, 
taking into accou nt the upper limits on the latter us ing the 
ASURV package llLavallev. Isobe. fc FeigelsonI Il99'^ . They 
also found a significant correlation between the specific SFR 
(sSFR) and X-ray luminosity, taking into account both up- 
per limits and a possi ble partial correlation with the redshift. 
iRosario et al.l (l2012li used the GOSMOS-GOODS- North 
and South X-ray selected sample. They study log{vLv) at 
60/rm vs. Lbol finding a significant correlation between Li, 
at 60gm and Lbol for moderate redshifts {z < 1) and high 
luminosities. 

Interestingly, all our objects are in the redshift range 
(1 < 2 < 3) studied by IPage et al.l ll2012h . They studied 
Herschel SPIRE observations of the CDF-N field (within 
the HerMES project) and found evidence for star formation 
(from 250 fim detections) in some X-ray detected AGN with 
10^® < Lx, 2 -io < 10“*^ erg/s, but lower star formation in 
the 10'*"^ < Lx, 2 -io < 10"*® erg/s luminosity bin. They took 
this as evidence for suppression of SFR in the most luminous 
AGN at that epoch in the Universe. Since all our objects are 
very significantly detected in 250 gm and have Lx, 2 - 10 > 
10"*^ erg/s (see Fig.|4]), in this scenario these highly luminous 
QSOs would not yet have managed to switch off SF. 

In this controversial context, we have revisited the 
\og{SFR) vs. log(I/x, 2 -io) correlation using both their sam- 
ple and a joint samp le with our sources (including data from 
llStevens et aLll^OOSl U. We have tested for a “hidden” corre¬ 
lation with redshift (specifica lly with the lumino s ity dis tance 
log(dL)) using the method in lAkritas fc Sieberd (Il996ll . who 
give their significance in terms of a ratio between the gen¬ 
eralized Kendall’s r and its dispersion a. 

We first started with just their sample, finding t/ a = 
5.1. We tested this significance against simulated samples 
of sources with mutually uncorrelated SFR and X-ray lumi¬ 
nosity values, but both correlated with redshift. Briefly, we 
found the constants K that best reproduced the log(y) = 
K -|-21og(di,) relations, with Y = log(SFi?) (using only de¬ 
tections) and Y = log(Lx, 2 -io), estimated the rms around 
these relations for several ranges of redshift, and then cre¬ 
ated 10000 samples of sources keeping the redshifts of the 
observed sources and simulating the X-ray luminosities and 
the SFR with the above “calibrations”. For upper limits in 
SFR we kept the observed upper limit but randomized the 
X-ray luminosity a s above. Thes e simu lated samples keep 
the statistics of the iRovilos et al.l (l2012l l sample but do not 
have any real correlation between SFR and X-ray luminosity. 
We found that 653 of those had r/cr > 5.1, so we conclude 
that the real significance of the correlation in their sample 
is about ~ 1 - 653/10000 = 93%< 2cr. 

We now wish to add our sample to this relation to 
check its infiuence. We nee d to use our full parent sam¬ 
ple in IStevens et al.l (l2005ll . including the non-detections 
at 850 fim. Note that there are three more “detections” 
in that paper compared to this paper: RX J1107-I-72 (be¬ 
cause it is radio-loud and hence its FIR emission originates 
in the AGN), RX J0943+16 and RX J1104+35 (these two 
are < 3cr significant). Th e 0.5-2 keV luminosities in that pa¬ 
per JSteveris et al.l[2005l i have been converted to 2-10 keV 
luminosities assuming Li, oc i'~^ (we have checked that this 
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method agreed w ell with the 2-10 keV luminosities from 
dPage et alJ l201lh used elsewhere in this work). For the 
five common sources, we have fitted the best multiplicative 
constant between t he 850 /rm-derived FIR luminosities in 
IStevens et al.l ll2005l ) and our SED-derived FIR luminosities, 
taking into account the uncertainties in both quantities. Us¬ 
ing this multiplicative constant we have then derived “cor- 
re cted” FIR luminosit ies and SFR (EQ. [8| for all sources 
in IStevens et al.l (|2005l i sample (green points in Fig[4|). The 
linked points show the magnitude of the re-scaling for the 
five common sources (red -this paper-, green -rescaled-). 
The rescaled full p arent sample has been joined to the 
iRovilos et al.l (120121 1 sample to study the X-ray luminosity- 
SFR correlation. Again, we have created 10000 random sam¬ 
ples as explained above (re-calibrating the SFR and X-ray 
luminosity correlations with redshift and their rms). The 
joint sample gave r/cr = 5.61. Only 23 out of 10000 ran¬ 
dom simulations showed higher values, so the significance 
is ~99.8%> 3 (t. This joint correlation is more significant 
probably because our five detected sources extend the X- 
ray luminosity range by almost an order of magnitude. 

At face value, it appears that there is a significant cor¬ 
relation between the growth of galaxies via star formation 
and the growth of their central SMBH via accretion in the 
joint sample. However, a number of caveats are in order to 
interpret the observed correlation. First, given the very dif¬ 
ferent selection functions of the samples involved (between 
them and among their constituent surveys at different wave¬ 
lengths) and the small numbers of sources involved, it is 
very difficult to assess the significance of the different re¬ 
sults in terms of the full population of black-hole-growing 
and/or star-forming galaxies at that epoch in the Universe. 
Large samples of objects at the relevant redshifts with well- 
controlled selection functions are needed to appraise this 
crucial issue. Furthermore, our SED-derived SFR cover very 
well the FIR rest-frame range even at our highest redshift (as 
they include observed frame points between 100 and 850 /rm) 
and the “corrected” ones used for the correlation come orig- 
inally from obser ved-frame 850 fim observations (similar to 
iLutz et aD 200 81. In contrast , the highest observed w ave¬ 
length in Royilos et al.l (l2012l 'l , iMullanev et al.l (|2012| ~) and 
iRosario et al.l 1 20121 1 is 160 fim, which only covers the short- 
wavelength side of the FIR bump at redshifts above 1, with 
the ensuing uncertainty in the FIR luminosities, despite the 
careful SED fit. 

Previous works have found some evidence of a corre¬ 
lation between the SF of the host and the AGN obscura¬ 
tion by neu t ral g a s in the X-ravs dA exander et al.l 2QQ^ 


Bajjer et al.l 12002 . iGeorgakakis et al.l 20041 . Rovilos et alJ 
200111 . Later studies with deeper surveys both in the X-ray 


and infrared ranges have failed to reprodu ce these results 
fe.g iRosario et al]|2012l . rRovilos et al.ll201^ 1. We have there¬ 
fore looked instead for a c orrelation betwee n the ionised col¬ 
umn density ('from l^ge et al]l201ll l and SFR in our 

sources (Fig. [5]), finding a tentative positive correlation be¬ 
tween these parameters. This is interesting, since it would 
imply a coupling of the ionized gas absorbing the X-rays 
at the scale of the accretion disk or the BLR with the gas 
forming stars in the host galaxy bulge, about three orders of 
magnitude farther away. At face value, th is would be com - 
patible with a positive feedback scenario, (iKing et al.ll201 j l 
in which the ionized outflowing gas would trigger star for- 



Figure 4. Star-formation rate versus X-ray luminosity for our 
five quasars (red solid diamonds, with values as reported in this 
paper). For comparison, w e also show X-ray selected AGN in 
deep Herschel surveys from iRnvilos et*^ I I2OI2I : 99 FIR-detected 
(empty dots) and 32 with only FIR upper limits (empty trian- 
gles). We have also included data from our full parent sample in 
IStevens et al.ll2005l. modifying the X-ray luminosities and SFR to 
a frame coherent with the data in this paper (see text): detected 
sources (green solid circles, joining with a segment the points cor¬ 
responding to the common sources) and upper limits (green solid 
triangles). 


1600 

1400 


1200 - I- "-1 

I_,,_I 

g 1000 ■ i 

s 

800 - I— ' —I 

cd 

600 - 

400 - I _ I 

22.0 22.5 23.0 23.5 
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sample. ffrom IPage et ahlDoilh . 


mation in the interstellar medium of the host galaxy, with 
the highest column density gas corresponding to stronger 
feedback. 

Alternatively, the AGN may also be ionising gas at kpc 
scales, co-located with the SF gas, so the above correlation 
would just be a consequence of the gas of higher density 
forming stars more intensely. Testing this intriguing possi¬ 
bility (with implications on positive and negative feedback) 
would need larger samples with a better determination of 
the location of the ionised gas. 
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4.4 Evolutionary status 

We now turn to the fate of the objects in our sample. As¬ 
suming constant accretion rates, we have estimated in Sec¬ 
tion EH their BH mass-doubling times r and the time they 
would take to reach the maximum BH mass observed locally 
Tmax- We note that RX J0057, RX J0941 and RX J1249 are 
already within a factor of 2-3 of this maximum local BH 
mass (~ 2 X 10^° A4q), so their BH growth phase is expected 
to finish within one to three mass-doubling times at most, 
i.e. 700, 1200 and 80 million years respectively. Given the 
expected li fetime of an active Q SQ phase of about 200 mil¬ 
lion years (|HoDkins et aLllioOSh . RX J1249 would have time 
to reach that maximum mass, while the other two objects 
would stop at a lower, but not much lower, mass. The BH 
in the centre of RX J1633, the least massive in our sam¬ 
ple (but already with a considerable mass of about 10^ Mp, 
large r than about 60% of the objects in iMarconi fc Hunu 
l2003lj . is growing quite fast, doubling its mass in 40 million 
years, so it could increase its mass by a factor of a few within 
the fiducial QSO lifetime above. RX J1218 has a intermedi¬ 
ate BH mass in our sample (about 10 times lower than the 
local maximum), but its accretion rate is the lowest in the 
sample, so its mass-doubling time is 3 Gy. Therefore, unless 
something happens to re-kindle the AGN at a later time, it 
is unlikely to grow much more. 

Without a determination of the masses of the host 
galaxies of our objects we cannot perform a similar exercise 
for the expected lifetime of the starburst phase. Unfortu¬ 
nately, the UV-to-MIR range is completely dominated by 
the AGN light. We are trying to secure mm observations to 
estimate the host galaxy masses directly. Nevertheless, we 
can argue that the SFRs observed are very high (albeit with 
considerable uncertainties), among the highest observed at 
the relevant redshifts, so they are unlikely to be maintained 
for a long time. 

Considering further episodes of strong galaxy growth 
with little AGN growth does not help to escape that the 
galaxies should be already formed, since either the host 
galaxy has to swallow whole fully-formed large-mass galaxies 
without forming stars (and with insignificant BH masses), 
or the ensuing SFR would have to be several times larger 
than the already extraordinary ones found in our objects 
(while at the same time avoiding a significant infall of gas 
to the galaxy nucleus to make sure to starve the massive BH 
already in place). 

In an independent line of argument, in the event that 
there is a positive feedback at work, as discussed in Section 
EH then the bulk of the gas mass in the host galaxies could 
be about to form stars pressured by the now maximal AGN, 
and so perhaps that could allow the BH-bulge relation to be 
reached quickly. However, this would be surprising, because 
the ~ 10^^ Mq of stars in an ellipti cal today typically seem 
to have formed at higher redshifts (iDaddi et al.ll2005ri . 

In the context of co-evolution of AGN and their host 

J alaxi es, e.g. the recent recently proposed by iLapi et ^ 
201411 . our objects, with 0.04 < Lfir/ Lbol < 0.81 would 
be in a stag e when the FIR- luminous phase is close to end 
(Figure 15 in lLapi et al.l20f3l . in qualitative agreement with 
our conclusion that their host galaxies are already mostly 
formed. 


5 CONCLUSIONS 

We have studied a sample of X- ray-obscured QSOs at 2 ~ 2 
with strong submm emission (IStevens et al.l 120101 1 . much 
higher than most X-ray-unobscured QSOs at similar red¬ 
shifts and luminosities, which, however, represent 85-90 per 
cent of the X-ray QSOs at the epoch. We have built X-ray- 
to-FIR SEDs for each object and we have fitted them us¬ 
ing different models (a direct AGN accretion disk, a torus- 
reprocessed component and a star formation component). 
We confirm that direct AGN, reprocessed AGN and SF com¬ 
ponents are needed to correctly characterize the SEDs of our 
objects. We have used these fits, together with our previous 
determinations of their X-ray luminosities, to estimate the 
total direct and reprocessed AGN luminosities and the lu¬ 
minosity associated to the star formation, as well as other 
derived physical quantities. 

We confirm the presence of strong FIR emission in 
these objects (well above that expected from plausible 
AGN emission models) which we attribute to SF at the 
ULIRG/HLIRG level with SFR~1000 Mq/j. Their associ¬ 
ated greybody temperature values are close to those of Sub¬ 
millimeter Galaxies (SMGs). They have dust masses around 
10’^ Mq. 

We have found a just over 3 -ct significant correlation be¬ 
tween the SFR and the X-ray lumino si ty wh en joining our 
sample with that from iRovilos et al.l (l2012l l. This is usu¬ 
ally taken as evidence for joint AGN and galaxy growth, 
but the differences between the techniques used to detect 
and characterize SF in those two samples detract from this 
otherwise exciting, interpretation. Our objects fall in the 
high luminosity end of large samples from deep surveys 
llRovilos et al.l (20121 : I Rosario et al.l (20121') , but have similar 
properti es to other samp les of high z, high FIR luminosity 
objects (iLutz et al.ll200^ . However, RXJ1249 stands out in 
all cases: it is one of the most luminous objects known (bolo- 
metric luminosity 10"*® L©) and it has an exceptionally 
high torus-to-bolometric luminosity ratio for its luminosity. 

Comparing their AGN reprocessed and direct emis¬ 
sion, we have obtained the ratios between the reprocessed 
emission from the torus and the total AGN bolometric lu¬ 
minosity ~ 0.3 — 0.9, higher than QSOs of similar lumi¬ 
nosities at 2 < 1.5. The Eddington ratios of our objects 
(0.1-0.6) are common for their redshift range, except for 
RX J1633 (our highest redshift object) which has a value 
of that ratio close to unity. Overall, the bolometric correc- 
tion s of our objects do not fit well with those of oth er stud¬ 
ies llVasudevan fc Eabian|[2009l : IMarconi et ^|2004| ~). show¬ 
ing higher bolometric luminosities compared to their X-ray 
luminosities. 

We have found a tentative positive correlation between 
NHion and SFR, perhaps indicative of positive feedback be¬ 
tween the X-ray and UV-detected ionized outflowing gas and 
the interstellar medium of their host galaxies. 

The black holes powering our QSOs are very massive at 
their epoch, ~ 10® — 10^°Mq (measured from broad emis¬ 
sion lines in their optical-UV spectra). We have calculated 
their mass-doubling timescale r and the time to reach the 
maximum BH mass observed locally, concluding that they 
can not grow much more. A further hint in this direction 
comes from the high Eddington ratio of RX J1633 which, 
according to iKellv fc ShenI ( 2013ll should persist only for a 
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very brief period of time. RX J1249 could become one of the 
most massive objects known. 

We do not know the masses of their host galaxies, but 
their black hole masses and their high SFR lead us to con¬ 
clude that they are already very massive or they would not 
have enough time to reach the local bulge-to-black-hole-mass 
ratio. This is also in agreement with recent models of AGN- 
host galaxy co-evolution. 

Direct determinations of the gas mass and of the mass of 
the host galaxies our QSOs are needed to have a better grasp 
of the nature and evolutionary status of these exceptional 
objects, and hence to understand their role in the disputed 
landscape of co-evolution of galaxies and AGN. 
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APPENDIX 


We now discuss briefly the fit results for each source: 

• RX J0057: We have cho sen the first torus model family 
from iRoseboom et al.l ll2013ll (11 members), since the other 
fits gave values at least 40 per cent higher. 

• RX J0941: We have s el ected the hrst torus model fam¬ 
ily from [Roseboo^^et_alJ_|201^ and the empirical tem¬ 
plate from Rowan-Robinson et al.l (120081') (22 m embers). The 
other torus models from Roseboom et al.l (120131 1 have val¬ 
ues at least 30 per cent higher. 

• RX J1218: Contrary to the other sources, the best-fit 
models were not visually good, since the y^ value was dom¬ 
inated by the hts to the optical-to-MIR region (where the 
error bars are smallest), leaving the SF bump badly repre¬ 
sented. For each diskJ-torus-l-SF combination we fixed the 
model parameters to their values best fitting the overall 
SED, and re-calculated the y^ value restricted to the rest- 
frame band 36-190 pm (5 points). We found that this re¬ 
stricted y^ had a clear best-fit peak with ^ 12, including 
20 disk-|-torus-|-SF combinations, with the rest of the com¬ 
binations extending in a tail towards higher restricted y^ 
values. We have chosen these 20 combinations as our best-fit 
“family”. Their torus co mponents include both the hrst and 
second torus models of IRoseb oom et"^ (l201 3li. as well a s 
the empirical template from iRowan- Robinson et al.l ll2008li . 

• RX J1249: This case is very similar to RX J0057: 

we have selected th e hrst torus model family from 
IRoseboom et ah] (l2013l i , for which the y^ values were at least 
40 per cent lower. _ 

• R X J1633: The third torus model from IRoseboom et al.l 
(120131) produced a well-grouped family (11 members) of 
best-hts with low y^ values. There were a few scattered 
best-hts with lower y^ using other models but they showed 
systematic residuals around rest-frame 30-60 pm, with SF 
components taking over the torus contribution in that band. 
We have chosen that family as our preferred ht. 

From the result s, see FigurelH above, it is clear that the 
hrst torus model of IRoseboom et ^ (l2013l f is preferred by 
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two of our sources (RX J0057 and RX J1249), and accept¬ 
able for another two (RX J0941 and RX J121 8), fo r which 
the empirical model of iR.owan-Robins on et af is also 

admitted. The second torus model of Roseboom et al.l (l2013l i 
only appears among the best-f its in one case (RX J121 8). Fi¬ 
nally, the third torus model in iR.oseboom et aP (l 201 . 8 l i is the 
best representation for RX J1633 

This paper has been typeset from a fQiiX/ IXQ 5 X file prepared 
by the author. 













